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Abstract. We present diagonal and off-diagonal elements of the permittivity 
and conductivity tensor up to the second order in magnetization for cubic crystals. 
We express all tensor elements as a function of a general sample orientation, 
for arbitrary magnetization direction and for (001), (Oil) and (111) surface 
orientations. Finally, we discuss ways, how to extract quadratic elements of the 
second order permittivity tensors for different sample surface orientations. 
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1. Introduction 

There is a vast number of physical phenomena, which scales with quadratic form of 
magnetization (also called effects quadratic in magnetization). Well-known examples 
are anisotropy magnetoresistance (AMR) [TJ [2] , quadratic magneto-optical Kerr effect 
(QMOKE) 0S1|5], longitudinal Hall effect j6], magneto-refractivity [3 [8], Cotton- 
Mouton effect (also called Voigt effect) or X-ray magnetic linear dichroism (XMLD) 
[51 [TO] (in some literature also called Hubert-Schafer effect [TO]). 

Although all those effects are not usually understand as a single phenomena, 
basically all of them originate from non-zero conductivity er^ (or permittivity ey) 
tensor, defines as = UijkiEjMkMi, where ji is the vector of current density, aijki is 
the fourth-rank tensor, being quadratic (second order) in magnetization (components 
Mj., Mi) and Ej is vector of the electric field. The difference between different 
effects mentioned above is than basically (i) photon energy range and (ii) whether 
the measured quantity is related with the diagonal or the off-diagonal part of the 
conductivity or permittivity tensor. Table [l] overviews techniques to study the off- 
diagonal permittivity or conductivity tensor in the first and in the second order in 
magnetization, and in the difference of the diagonal terms. 

With some generalization, basically every magneto-transport phenomena is 
expected to have its second-order counterpart. Hence, existence of effects such as 
quadratic spin-Hall effect or quadratic magneto-Seebeck effect can be envisaged, 
although those effects have not been demonstrated experimentally neither by ab-initio 
calculations yet. 

Within this Article, based on symmetry arguments, we determine the diagonal 
and the off-diagonal elements of the permittivity tensor in cubic crystal up to 
the second order in magnetization. We express all tensor elements for arbitrary 
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conductivity 


magneto-optics 


X-ray 


off-diagonal elements lin- 
ear in M, s\f, i =/= j 


Hall effect 


MOKE 


XMCD 


off-diagonal elements 
quadratic in M s\j , i J 1 j 


longitudinal Hall 
(quadratic Hall*) 


QMOKE 


quadratic XMCD* 


difference of diagonal 
terms, e£ 2) - ef) (i ^ j) 
being quadratic in M 


AMR 


Voigt effect 
(Cotton-Mouton effect) 


XMLD 



Table 1. Techniques used to study first and second order effects in magnetization. 
Techniques denoted by * are not experimentally demonstrated yet. 



magnetization direction and for (001), (Oil) and (111) surface orientations of the 
sample. We also distinguish two different kinds of symmetries in cubic crystals, (i) 
the crystals with the symmetry classes 23 and m3 and (ii) after adding the symmetry 
operation the rotation about 180° about xy axis, the cubic crystals with the symmetry 
classes 43m, 432 and m3m are described. 

2. Motivation 

Within the complex representation of the Maxwell equations, the complex (relative) 
permittivity tensor and the complex conductivity tensor tr^ describe the same 
physical phenomena, with conversion between them being (in SI) £y = + i/Jij/u, 
where 8™ is Kronecker delta and u is the photon frequency (here, ay is expressed in 
units [s -1 ]). Hence, from symmetry point of view, the symmetry form of and <7jj 
are exactly the same. Therefore, in following we express tensors' symmetry form solely 
for permittivity. The form of conductivity tensor is obtained simply by interchanging 
variables o~%j O As we express the tensors' symmetry form in permittivity e^ , 
we later call this tensor by its historical name as magneto-optical permittivity tensor, 
although, as told above, the presented symmetry analysis are not limited to magneto- 
optical investigations. 

Those symmetry arguments analysis are important in following cases: 

• Experimental evidence of the second order effects may be difficult, as those 
effects can be small or they can be mimic by artefacts of the experiments. 
Therefore, to confirm that the measured effect is a quadratic effect, one should 
measure interplay of magnetic field direction, crystallography axis direction and 
experimental geometry. For example, when one founds geometry where first order 
effects vanishes and when the behaviour of measured effects agrees with the form 
predicted by the symmetry arguments, then the presence of quadratic effects is 
clearly proven [H[S]. Therefore, detailed understanding of symmetry of quadratic 
effects helps to optimize measurement geometry for each experimental technique. 

• Furthermore, those combined measurements allow to determine quadratic 
elements of the second order permittivity tensors (in following denoted as G44, 
AG, Ar). The symmetry arguments are important to determine experimental 
geometry, in which those elements can be extracted. This opens new spectroscopy 
fields where spectra of G44, AG, Ar quadratic elements are determined. Up to 
now, no spectra of quadratic elements are available and the complex values of 
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G44, AG were determined only for Fe(OOl) for a single wavelength [llj . 

• The symmetry arguments presented here are also important for ab-initio 
calculations, to verify correctness of the calculated permittivity or conductivity 
on applied direction of the magnetization as well as to express ab-initio spectra 
of G44, AG, Ar quadratic elements. Nowadays ab-initio codes calculate total 
permittivity or conductivity, as a function of the applied direction of the magnetic 
field. Hence, in order to separate the quadratic elements from the ab-initio- 
calculated total permittivity or conductivity, symmetry arguments predictions of 
the dependence of the tensors on magnetization direction must be used and fitted 
to the ab-initio values of the tensor elements. 

• Last, note that definition of second order permittivity elements and their 
dependence on magnetic field and sample orientation contain several sign 
conventions. However, to correctly compare extracted second order permittivity 
tensors from both experiment and ab-initio calculations, the sign conventions 
must be handled carefully and correctly. The tensor transformation presented 
within this Article were done with a careful definition of all sign conventions. 
Therefore, this Article may also establish the sign conventions for second order 
effects. This will be particularly important in expected spectroscopies of the 
second order permittivity tensors. 



3. Magneto-optical permittivity tensor 

In magneto-optically active materials we can evolute permittivity tensor with respect 
to the components of magnetization M acting on the material [3] 



-(0) 



deij_ 
dM k 



4 0) + K ijk M k 



Mk + l 

M=0 Z 

GijkiMkMt 



d 2 e t 



dM k DM, 



M k Mi + 



(1) 



M=0 



where — 



£ji are components of permittivity tensor of material without 



acting magnetization M. Furthermore, e 



(1) 



KakMh denotes contribution to 



the permittivity tensor linear in magnetization M where a third rank tensor 

is called linear magneto -optical tensor. Contribution to the permittivity tensor 
(2) 

e\.' = GijkiMfrMi is quadratic in magnetization M and is described by the fourth 
rank quadratic magneto-optical tensor Gijki- 

Magneto-optical permittivity tensor must fulfil Onsager relation 



e tf (M) = e 3 -i(-M) 



(2) 



(21 

Hence, the quadratic magneto-optical contribution to the permittivity tensor must 



be symmetric [5] 



£ (2) 
Gijki 



-(2) 
"ji ' 



which implies 



(3) 



Gjiki — Gijik — Gjuk 
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Thus, the quadratic MO permittivity tensor can be rewritten in simplified form [3j 
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(4) 



being a general expression of the second order MO permittivity tensor. 
3.1. Quadratic MO tensor in cubic crystals 

In our work we are focusing on cubic crystals. Let's suppose that coordinate system 
corresponding to the principal axes of the crystal (Fig.[l]). Then we can write quadratic 
MO tensor in cubic crystals in general form (the crystal classes 23 and m3) [T2"] 
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Figure 1. Principal axes of cubic crystal. 
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(5) 



where short index notation of Gj,-fc; was used. Adding the symmetry operation the 
rotation about 180° about xy axis we have for cubic crystals with the crystal classes 
43to, 432 and m3m [H] 
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(6) 
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x=001 



Figure 2. Transformation of vector's components under the orthogonal 
transformation of coordination system. 



Adding further symmetry operation then describes isotropic medium [12] 
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(7) 



For interpretation of further calculations, it is convenient to introduce new terms AG 
and Ar 

Gl2 + G21 



AG — Gn — 2G44 — 
_ G12 - G 2 i 



(8) 



Using those terms, the components G12 and G21 are then expressed as G12 = 
Gn - 2G 44 - AG + Ar and G 2 i = Gn - 2G 44 - AG - Ar. In cubic crystals with 
higher symmetry (the crystal classes 43m, 432 and m3m) the parameter Ar vanishes, 
Ar = 0. Furthermore, in case of isotropic medium also the parameter AG vanishes, 
AG = 0. 



4. Transformation of vectors and second rank tensors under orthogonal 
transformation of coordination system 

Now we want to describe how the quadratic MO permittivity tensor is changed under 
the rotation of coordination system. The original coordination system (later called 
'crystal coordinate system') defines xyz axes being parallel to the principal axes of 
crystal. 

At first let's suppose that we transform crystal coordination system according to 
an orthogonal transformation given by matrix R (see Fig. [2| . Having a vector A in 
crystal coordinates xyz, then the same vector in new coordinates x'y' z' is expressed 
by components A' [13 (in Einstein notation) 

A[ = RijAj (9) 

or written equivalently as A' = R • A in matrix notation. Next, we transform the 
components of a second rank tensor [13] 



(10) 
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or equivalently written as e' = R • e ■ R . in the matrix notation. Now we perform two 
orthogonal transformations of coordination system in succession. First we transform 
coordination system according to the matrix R*- 1 ' and then according to the matrix 
R( 2 ). Then the coordinates of the vector A in new coordinate system are given by 
complete transformation Rij = R^} R& (i.e. in matrix form R = R( 2 )RW). 



5. Crystal in a hypothetical experimental setup 

Now, let us place the cubic crystal (i.e. the sample) into a hypothetical setup, described 
by cartesian axes x"y" z" . Then, the parameters related with the setup, such as e.g. 
the external field direction, plane of incidence of light, flow of charge current, sample 
surface etc. are described within this setup cartesian system. Furthermore, without 
loss of generality, let us attribute z"-axis to be normal to the sample surface and hence 
x" and y" directions are within sample plane. 

As wrote above, the main crystal axis (100), (010) and (001) are related with 
xyz coordinate system called as "crystal coordinate system" , in which the quadratic 
response is described by Eqs. ([5]-[7])- 

In following Sections, we describe quadratic response of cubic crystal, where 
sample surface (i.e. sample normal) has orientations (001), (011) or (111). Hence, 
as the sample surface orientation is fixed within crystal coordinate system, the only 
free parameter is rotation of the sample around the sample normal, i.e. rotation around 
z" axis by an angle 7. In real experimental investigations, it corresponds to rotation 
of the sample in case of optical investigations, or to change charge current direction 
in case of transport investigations. The positive 7 direction is defined as follow: when 
looking into the positive z" direction, the rotation is counter-clock- wise (Fig. [3]) . 

Finally, recall that the relation between both coordinate systems is given by 
matrix R as in Eq. ^. 

6. Quadratic MO permittivity tensor for cubic crystal with (001) 
oriented surface 

Let's suppose a cubic crystal grown with (001) surface direction. Hence, the setup 
coordinate system and the crystal coordinate system have equal z-axes, z = z" 
pointing in sample surface normal direction. Then, the rotation around angle 7 is 
described by transformation matrix R 2 (7) (Fig. [3]) 



R(°°D=R,( 7 ) 



cos 7 sin 7 
— sin 7 cos 7 
1 



(11) 



Before expressing the elements of permittivity tensor given by (|5jj7j) in the 

setup coordinate system x"y"z", £^ 01 \ first we need to express components of 
magnetization vector in the crystal coordinate system. It originates from the fact 
that during hypothetical measurements, we describe magnetization vector M" in 
setup coordination system x"y"z" ', but elements M k , Mi coming in Eqs. ^fff§ are 
components of the same magnetization vector M expressed in the crystal coordination 
system xyz. Recall, the relation of the magnetization vector between both coordinate 
system is (Eq. Q) 

M k =R^M' s '=R^M'J. (12) 
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z=001=z' 




x=100 



Figure 3. Transformation of coordination system for sample with (001) oriented 
surface. Presented angle 7 is positive. 



Hence the permittivity tensor in crystal coordinate system xyz with magnetizations 
from setup coordinate system x"y" z" is 



4 2 > = G^i^C^ 001 ^ 1 ^ 



(13) 



(2) 

Finally, transformation of the permittivity elements e\, into the sample coordinate 



system is (Eq. [10 



(001) _ R (001) p(001) (2) 
t ij — n ik 11 jl t kl 



(14) 



which leads to the final expression of the quadratic MO tensor in the setup coordinate 



(001) 



with magnetization direction also described in the setup coordinate 



system, e 
system 

£ (ooi) = ^^G^ooi)^^)-^ . (15) 

where recall Gumn are elements of the quadratic MO tensor, as described by Eqs.fJHJ- 
m). In Tabs. [2] and [3] are expressed calculated off-diagonal and diagonal elements 
of j as well as AMR-like parameter e^® — £t/y° 1 ' ) • In corresponding rows are 
separated terms proportional to elements of magnetization MyM'J, M"M'J , M"My 

and M" 2 y - M" 2 Z , M"\ ~ M" 2 X , M'\ - M'\. 

For comparison of experimental or ab-initio results with symmetry arguments, it 
is often useful to express magnetization vector M" in spherical coordinates 



M" 



M'i 



sin t) cos ip 
sin 9 sin ip 
cos 9 



(16) 



Then, all permittivity elements can be expressed to depend on angles p and 9 

as presented in Tab. [4] 
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7. Quadratic MO permittivity tensor for cubic crystal with (Oil) 
oriented surface 

Here, we discuss a cubic crystal grown with (Oil) surface orientation, i.e. the normal of 
the sample surface (corresponding to z' direction) is parallel to the [Oil] direction. To 
obtain new coordinate system x'y' z' we rotate old coordinate system xyz about angle 
— 7r/4 about x-axis (Fig. PO). This is described by transformation matrix R x (— 7r/4) 



R^-tt/4) = 



1 

y/2/2 -V2/2 
s/2/2 V2/2 



(17) 



z=001 




Figure 4. Transformation of coordination system for sample with (Oil) oriented 
surface. 



Furthermore, let us discuss rotation of the sample by an angle 7 around its sample 
normal, i.e. let us rotate x'y 1 plane around z' = z"-axis by the angle 7, providing setup 
cartesian system x"y"z" . This transformation from coordinates x'y' z' to coordinates 



x"y"z" is given again by R z (7) matrix (111. The total transformation matrix from 



xyz crystal coordinates to x"y"z" setup coordinates is 

R(° n > = R,( 7 )R B (-7r/4) , (18) 
so we can transform elements of quadratic MO permittivity tensor similar as in 



previous Section (Eq. (15)) 



£ i 3 - K ik H jl ^klmniO ms M s H nt M t ■ U 9 ) 

Similarly as in previous section, in Tab. [5] are calculated off-diagonal elements of 
s^j 11 ^ and AMR-like parameter eS; 11 ^ — £yy^ and Tab. ji] provides calculated diagonal 
elements of e^ 11 ^, in both cases as a function of quadratic form of magnetization 

components M'^M". Finally, Tab. [7] provides all elements of Ey 11 ^ expressed as a 
function of spherical coordinates of magnetization M", angles ip and 9, as defined in 
Eq. m. 
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8. Quadratic MO permittivity tensor for cubic crystal with (111) 
oriented surface 

Let us discuss the cubic crystal with (111) oriented surface. Then, the z'-axis, being 
perpendicular to the sample surface, is parallel to the [111] direction. To obtain new 
coordinate system x'y' z' we rotate crystal coordinate system xyz about angle a around 

i F) 

— zy-axis, where cos a = ^ and sin a = 

To calculate transformation matrix we use well-known transformation matrix for 
rotation about arbitrary angle a about arbitrary axis given by normalized vector 
n = (ni,n 2 ,n 3 ) 



z=001 z'=111 




x' 



Figure 5. Transformation of coordination system for sample with (111) oriented 
surface. 



R{a) 



cos a + n\(\ — cos a) nin 2 (l — cosa) + n% sina n\n^{\ — cos a) — n-i sin a 
^1^.2(1 — cosa) — n3sinct cos a + n|(l — cos a) 722713(1 — cosa) + m sina 
^1^-3(1 — cos a) + n2 sin a 712713(1 — cos a) — n\ sina cos a + rz§ (1 — cos a) 

(20) 

In our case (m, 712, 713) = (—^75, 0) and sina = i we have 

iii iii 1 



R(a) 



2 2V3 

71 



2^3 

j L_ 

2^3 

71 



t 

Y 3 



(21) 



Now we can rotate sample by an angle 7 around z' — z"-axis about arbitrary angle 
7 given by transformation matrix R z (7) (111. Then, the total transformation matrix 



from xyz crystal coordinate system to x"y" z" setup coordinate system is 

R< U1 > = R,(7)R(a) (22) 
and we can transform elements of quadratic MO permittivity tensor similarly as in 



Eq. (15) 



Jin) _ d(hi) R (iii)^ R( m ) _1 a r" p^ 11 ) -1 M" 

£,'-,' — Mil. Lrklmn-tl ms M s K nt M t 



-(HI 



(23) 

Similarly as in two previous cases, in Tab.[8]are calculated off-diagonal elements of 



and AMR-like parameter e^ -e^ and in the Tab.g are calculated diagonal 



elements of e 



(111) 



In the Tab. 



10 



are all elements of ejj expressed as a function of 



the spherical coordinates of magnetization vector M" (16). 
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9. Discussion 



To determine quadratic elements, Gy, in case of optical investigations, one can vary 
field direction (described by angles tp and 9) or one can rotate sample around its axis 
(described by an angle 7). In case of the conductivity measurements, the equivalent of 
the sample rotation is in-plane rotation of the in-plane current with respect to crystal 
axis. Namely, it is equivalent to the case when current flows in a fixed direction with 
respect to the setup cartesian system, and only crystal itself is rotated. Therefore, from 
symmetry point of view, following discussion is valid for both optical and transport 
investigations. 

Now, let us discuss possibility to separate Gij quadratic elements from 
measurements or ab-initio calculations. We will discuss hypothetical measurements or 
calculations for three magnetization scans. First scan, 'in-plane magnetization scan', 
keeps magnetization within a;"i/"-plane, i.e. the magnetization orientation is described 
as 8 = Tt/2 and ip varies. This scan is particularly easy to obtain by experimental 
means, as well as (for 7 = 0) by ab-initio calculations. Second magnetization 
scan, 'out-of-plane-I magnetization scan' goes from ^''-direction (i.e. from out-of-planc 
magnetization, i.e. from 8 = 0) through (111)" direction towards (110)". This scan 
corresponds to tp = 7r/4 with varying 8. This scan is particularly is important for 
ab-initio calculations as in case of (001) oriented crystal with 7 = 0, as this scan goes 
through all three types of high symmetry points in cubic crystal, namely (001), (111) 
and (110) [13]. The last magnetization scan to be discuss is called 'out-of-plane-II 
magnetization scan' being defined as tp = and 9 varies, i.e. it goes from out-of-plane 
magnetization orientation (001)" to in-plane orientation (100)". 



9.1. Magnetization scans for (001) oriented surface 

First we treats in detail (001) oriented sample. The dependence of all permittivity 



elements for each magnetization scan is presented in Tab. 11 (i) It reveals that the 
term Gn is not accessible by any means, as it stays constant for any sample and 
magnetic field orientation and hence it can not be extracted, (ii) The in-plane scan 
can provide values of all G44, AG and Ar from both the off-diagonal element (Hall- 
like) investigations and diagonal element (AMR-like) investigations. The extraction 
is based on different periodicity of the quadratic elements on 7. However, in case 
of Ar, this may be difficult to extract Ar ~ sin 27, as experimental artefacts (such 
as misaligned sample) may have the same symmetry, (iii) The most reliable way for 

extraction of Ar is from e£ 01 \ as it varies solely on Ar through parameters tp and 7. 
However, from experimental point of view, this element is very difficult to obtain by 
optical investigations and in case of the conductivity investigations, it requires current- 
perpendicular-to-plane geometry, (iv) In case of out-of-plane magnetization (8 = 0), 
the off-diagonal element eiy 01 ^ simplifies to ei°° 1 ' ) = Arsin27. This provides probably 
the simplest way to determine Ar by optical investigations, simply by observing 
anisotropy of the polar Kerr effect (or Hall effect or XMCD) on sample orientation 7. 



(v) Finally, Table 11 demonstrates, that in case of out-of-plane scans, the periodicity of 
different permittivity elements on 7 is similar as in case of in-plane magnetization scan. 
Also, no simple separation of quadratic term appears. Therefore, from experimental 
point-of-view there is no particular advantage in quadratic parameter separation to 
apply scans varying magnetization between out-of-plane and in-plane directions. 
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9.2. In-plane magnetization scan for (Oil) and (111) oriented surface 

In previous paragraph, we have demonstrated that there is no particular advantage 
from experimental point of view to use out-of-plane magnetization scans to separate 
contributions from different quadratic elements. Therefore, in case of (Oil) and (111) 
oriented surfaces, we discuss dependence of permittivity elements only for in-plane 
magnetization scan, presented in Tabs. p~2] and [13} respectively. 



In case of (Oil) surface probed with the in-plane magnetization scan (Tab. |l2|), 



one can see that (i) both Hall- like term (e^y 1 ^) and AMR-like term (el™ 1 '' — £yy) 
depend solely on G44 and AG and their value can be extracted by both varying 
magnetic field ip or sample orientation 7. Like in (001) surface oriented crystal, G44 
contribution is independent on sample orientation 7 whereas AG contribution varies 
with 7. Interestingly, in later case the dependence of AG has periodicity on both 27 
and 47. (ii) Ar is accessible solely through Sy ^ and £m terms. However, when 
magnetization goes out-of-plane, dependence of all terms on Ar appear (see Tab.[9|. 

In case of (111) oriented surface probed with in-plane magnetization scan 
(Tab. |13[ ), remarkable point is that the periodicity is either 2ip (when rotating 

magnetization) or 37 (when rotating sample). For eiy 11 ' (Hall-like), ei^ 11 " 1 — e^y 1 ^ 
(AMR-like) permittivity elements and for all diagonal permittivity elements, all three 
quadratic elements G44, AG and Ar contributions are either constants or they 
depend on magnetization orientation as 2tp, all of them being independent on 7. 
Remaining permittivity elements £xz and s^z^ depends only on AG and Ar and 
both contributions in a similar dependence on 7 and ip. Therefore, all it makes 
separation of values of quadratic elements from (111) oriented crystal difficult. 

10. Conclusion 

Using symmetry arguments, we determined diagonal and off-diagonal elements of the 
permittivity tensor up to the second order in magnetization for cubic crystals. We 
express all tensor elements for arbitrary magnetization direction and for (001), (011) 
and (111) surface orientations. Finally, we discuss a way to separate various quadratic 
elements of the second order permittivity tensors for different magnetization direction 
and different surface orientations in case of both optical and transport investigations. 
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Appendix A. Linear-in-magnetization permittivity and conductivity 
tensor for cubic crystals 

Although this Article treats permittivity/conductivity tensor in the second order in 
magnetization, let us shortly overview their dependence on the zeroth and the first 
order in magnetization. 
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M x = sin 8 cos <p 


Mj, = sin 9 sin ^ 


M z = cos e 
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Table Al. The first-order-in-magnetization permittivity elements eS.' as a 

function of magnetic field orientation for cubic crystal. ey' does not depend 
on crystal orientation. 



In case of the zero-order permittivity ef^ in the cubic crystal, the permittivity is 
independent on both crystal and field orientation, = e^Sij, where 5ij is Kronecker 
delta, and is scalar zero-order permittivity. 

In case of permittivity elements in the first order in magnetization, their 
dependence is summarized in Tab. Al written in similar notation as all other Tables 
introducing the second-order permittivity. Note that does not depend on crystal 



orientation, so Tab. |A1| is valid for any surface orientation of the crystal and any 
sample orientation 7. 
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